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The bulk photovoltaic effect is a long-known but poorly understood phenomenon. Recently, 
however, the multiferroic bismuth ferrite has been observed to produce strong photovoltaic response 
to visible light, suggesting that the effect has been underexploited as well. Here we present three 
polar oxides in the LiNbOa structure that we predict to have band gaps in the 1-2 eV range and very 
high bulk photovoltaic response: PbNiOs, Mgi/ 2 Zni/ 2 Pb 03 , and LiBiOs. All three have band gaps 
determined by cations with electronic configurations, leading to conduction bands composed of 
cation s-orbitals and O p-orbitals. This both dramatically lowers the band gap and increases the bulk 
photovoltaic response by as much as an order of magnitude over previous materials, demonstrating 
the potential for high-performing bulk photovoltaics. 


Photovoltaic effects have long been observed in bulk 
polar materials, especially ferroelectrics [H-Q- Known 
as the bulk photovoltaic effect (BPVE), it appeared to 
derive from inversion symmetry breaking. Despite in¬ 
tense initial interest, early explorations revealed low en¬ 
ergy conversion efficiency, in part due to the high band 
gaps of most known ferroelectrics. Additionally, despite 
several proposed mechanisms, the physical origin of the 
BPVE was unclear iil. 

However, recent emphasis on alternative energy tech¬ 
nologies and the observation of the effect in novel semi¬ 
conducting ferroelectrics (with band-gaps in the visible 
range) has renewed interest (Ml- Several studies have 
attempted to elucidate the various contributions to the 
photovoltaic response - bulk or otherwise - in ferro¬ 
electrics iB“ 23|. In particular, bismuth ferrite (BFO) has 
been found to generate significant bulk photocurrents; 
combined with its unusually low band gap of 2.7 eV, it 
has attracted a great deal of attention for its potential in 
photovoltaic applications [U [25 - l3lj| . The understand¬ 
ing of the fundamental physics behind the effect has ad¬ 
vanced as well; recently we demonstrated that the BPVE 
can be attributed to “shift currents”, and that the bulk 
photocurrents may be calculated from first-principles. 
The ab initio calculation of the shift current and sub¬ 
sequent analysis yielded several chemical and structural 
criteria for optimizing the response. These criteria have 
been used previously to modify or identify existing mate¬ 
rials with enhanced response 3, 3 - 341. In this work, we 
use these insights to propose several candidate bulk pho¬ 
tovoltaics with calculated response as much as an order 
of magnitude higher than well-known ferroelectrics, while 
having band gaps in or slightly below the visible spec¬ 
trum. Our results demonstrate that bulk photovoltaic 
response can be much stronger than previously observed, 
supporting the possibility of materials suitable for appli¬ 
cation. 


There are two figures of merit for evaluating the BPVE 
in a material: the current density response to a spatially 
uniform electric field, and the Glass coefficient [3|. The 
current density response is given by the tensor 

Jq{u}) = tTrsg(w)E°(w)A°(a;) 

O'r-sq(w) = e ^ J dklrs{n', n" , k; uj)TZq{n', n", k) 

n' ,n" 

where E° is the vector of the illumination field, and n' 
and n" index bands. Letting / denote filling, y the Berry 
connection, and (j) the phase of the transition dipole, the 
expression 
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describes the intensity of transitions, and 
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is the expression for the “shift vector”, which describes a 
distance associated with the excited carrier [1, 35|, and 
depends on the differences in the wavefunction centers up 
to a unit cell, as provided by the Berry connections, and 
the average separation in unit cells given by the transi¬ 
tion dipole phase derivative. Roughly speaking, the two 
terms X and TZ can be thought of as giving the number of 
carriers excited and the velocity of those carriers. We em¬ 
phasize that this mechanism is profoundly different from 
other photovoltaic effects; rather than relying on excita¬ 
tion of carriers which are then separated by an electric 
field, the carriers are electron/hole pairs in coherent ex¬ 
cited states that possess intrinsic momentum of opposite 










2 


sign. Crucially, this allows for arbitrarily high photovolt¬ 
ages; the Schockley-Queisser limit does not apply, a ma¬ 
jor advantage of BPVE. In particular, the open-circuit 
voltage is determined by the competition between the 
photocurrent and countervailing voltage-driven leakage 
that depends on the overall resistance of the sample [S^ ■ 
This sample dependence prevents straightforward calcu¬ 
lation. 

Determining the total current in a sample is compli¬ 
cated by the attenuation of incident illumination as it 
travels through the material. In the limit of a thick sam¬ 
ple that will completely absorb the illumination, the total 
current can be obtained from the Glass coefficient G 

\E°{uj)fyV = Grrg{oj)Ir{u;)W (3) 

^rr{u) 

where a is the absorption coefficient, and W is the sam¬ 
ple width. Thus, the current density tensor and Glass 
coefficient describe the response in the regimes of near¬ 
zero and near-infinite thickness, respectively. In practice, 
“infinite thickness” is on the order of microns, and total 
photocurrent is usually best described by the Glass coef¬ 
ficient. 

However, the Glass coefficient provides additional in¬ 
formation about the response. In the limit where e* ^ e’', 
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The Glass coefficient is therefore closely related to the 
weighted average shift vector, allowing us to estimate the 
contribution of both terms in the shift current expression. 
Shift current response was calculated as in Ref. [37j, 


from wavefunctions generated using density functional 
theory, with the generalized gradient approximation 
(GGA) and optimized, norm conserving pseudopoten¬ 
tials [H, Is^. The presented results exclude spin-orbit 
effects; calculations with and without spin-orbit were per¬ 
formed for both LiBiOa Mgi/ 2 Zni/ 2 Pb 03 and were not 
found to substantially influence the results. For BFO, a 
Hubbard U of 5 eV was used for Fe Sd, as in Ref. (ioj l. 
For PbNiO^, a Hubbard U of 4.6 eV was used for Ni 3d, 
as in Ref. [ill. QUANTUM ESPRESSO [H was used 


for the electronic structure calculations, and OPIUM was 
used to generate pseudopotentials. The Heyd-Scuseria- 
Ernzerhof (HSE) hybrid functional [i^ was used to com¬ 
pute band gaps, as it is known to frequently produce 
significantly more accurate values than GGA. These cal¬ 
culations were performed on 8 x 8 x 8 k-point grids, with 
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FIG. 1. The LiNbOa (LNO) primitive cell, overlaid with the 
pseudo-cubic perovskite cell. The direction of the polar dis¬ 
tortion is indicated by the black arrow. 


4x4x4 grids for the exact exchange HSE calculations. 
Band structures and density of states plots are gener¬ 
ated from GGA calculations, and reported HSE band 
gaps are the direct gaps. The present results are for the 
experimental structure in the cases of LNO, BFO, and 
PbNiOa, and computationally relaxed structures for the 
other materials. Structural relaxations and calculations 
of the shift current were performed at the level of LDA 
and found to vary minimally from the GGA results; due 
to the high expense of exact exchange calculations, the 
dense k-point grids required to converge shift current cal¬ 
culations cannot prese ntly be obtained using HSE, and 
scissor corrections 3 to the HSE gaps were applied 
to account for the dependence of Glass coefficient on fre¬ 
quency. 

Previously, we revealed the dependence of shift vec¬ 
tor magnitude on the chemical and structural properties 
of materials. Large shift vectors were characterized by 
valence and/or conduction states that are both stro ngly 
asymmetric and delocalized in the current direction (37j | . 
In this regard, many distorted perovskite {ABO 3 ) ferro- 
electrics are crippled by the presence of cations en¬ 
closed in octahedral oxygen cages. The conduction band 
edge is dominated by t 2 g-like d states that are largely 
nonbonding. Goupled with the tendency for d states to 
localize, the result is that both shift vectors and tran¬ 
sition response are very weak near the band gap. The 
delocalized Cg states are much higher in energy, effec¬ 
tively raising the energy threshold for significant BPVE. 


To overcome the weak BPVE response of oxides, 
we investigated systems that involve both large distor¬ 
tions to oxygen cages, (increasing the bonding character 
of any c?° states) as well as cations with less localized 
s and/or p states near the band edge [S^- It has already 
been noted that cations can dramatically improve the 
activity of photocatalysts [H3|- We found polar oxides 
taking the LiNbOa structure to be promising candidates. 
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PbNiOs 

Mgi/2Zni/2Pb03 

LiBiOs 

a 


5.63 A 

5.77 A 

5.67 A 

a 


57° 

57° 

56° 

A 

(2a) 

(0.0, 0.0, 0.0) 

(0.0, 0.0, 0.0) 

(0.0, 0.0, 0.0) 

B 

(2a) 

(0.214, 0.214, 0.214) 

(0.216, 0.216, 0.216) 

(0.213, 0.213, 0.213) 

0 

(6b) 

(0.830, 0.098, 0.415) 

(0.794, 0.128, 0.390) 

(0.798, 0.122, 0.405) 

P 


99 pG/cnP 

83 fiC/cm^ 

50 fiG/cnP 


TABLE I. The structural data for the three compounds presented here. PbNiOs and LiBiOs are in space group R3c, while 
Mgi/ 2 Zni/ 2 Pb 03 is in R3. However, the deviations of the coordinates from the R3c positions are miniscule (< 0.2%), so they 
are presented as such with Mg and Zn each occupying one site of the A position. Polarizations were determined based on a 
non-polar structure featuring the A-site atom coplanar with oxygen, and the B-site midway between oxygen planes 



FIG. 2. The current density response for BFO (GGA+U) is 


shown in (a) The Glass coefficient of BFO appears in (b) 


Only the response in the direction of material polarization 
is shown, for both perpendicular {xxZ) and parallel {zzZ) 
light polarization. Dashed lines appear at benchmark values 
of current density and Glass coefficient chosen to represent 
the maximum response of these materials. 


with cations and Bi®"*". This structure can 

also be obtained by distorting the perovskite structure 
rhombohedrally, and allowing polar distortions along and 
oxygen-cage rotations about (111). Notable ferroelectrics 
with this structure (but with cP cations) include LiNbOa 
(LNO) and BiFeOa (BFO). LNO is known for its large 
nonlinear optical response, and, often doped with iron, 
it was one of the first materials in which the bulk photo¬ 
voltaic effect was observed and studied [1, 51) 13 ■ How¬ 
ever, its bulk band gap is well outside the visible spec¬ 
trum [HsI . BFO has garnered much attention recently for 
its multiferroic behavior 5^ and low band gap of about 
2.74 eV [ 5 ^, which has led to explorations of its photo¬ 
voltaic response 21 - 24, . We have used BFO as a 

benchmars for the present study; as with the archetypal 
ferroelectrics BaTiOa and PbTiOs, its LUMO is dom¬ 
inated by cation d-states and yields a very similar re¬ 
sponse magnitude. 

We consider only current response in the direction of 
material polarization for both perpendicular {xxZ) and 
parallel (zzZ) light polarization, as these are the only 
tensor elements that can contribute to the response to 
unpolarized light. For ease of comparison, we mark base¬ 
line values reflecting the maximum response of our bench¬ 
mark, shown in Fig. [51 with a dashed, red line. These are, 
for the current density and Glass coefhcient, respectively, 


5 X 10 ^(A/m^)/(W/m^) and 5 x 10 ®cm/V. 


We have studied three materials taking the LNO struc¬ 
ture (Fig[T|): PbNiOs, Mgi /-^Zni /^PbO.^, and LiBiOs. 
The first has been synthesized |58l |. and the latter two 
are similar in composition to known materials. The struc¬ 
tural parameters and bulk polarizations are given in Ta¬ 
ble jT] The distortion from cubic perovskite is sufficiently 
strong that assignment of A- and B-sites is ambiguous; 
we have followed the assignment of Ref. for PbNiOs, 
but note that treating Ni as the A-site (reversing the 
orientation), the Wyckoff position of oxygen becomes 
(0.800,0.116,0.384), which is slightly closer to the cor¬ 
responding crystal coordinate of the other two materials. 
All three satisfy our requirements of low band gap, 
cations, and large polar distortions. Furthermore, as seen 
in Fig. |3l all three have qualitatively similar band struc¬ 
tures, featuring highly dispersive conduction bands, in 
contrast to the usual case of dP perovskite derivatives. 
As we will show, this arises due to unfilled s-like - rather 
than d-like - states composing the conduction band, and 
has profound consequences for the bulk photovoltaic re¬ 
sponse. 


PbNiOs has recently been synthesized [58| and ex¬ 
plored theoretically [^, [^. Like BFO, it is antiferro¬ 
magnetic with weak spin-canting, and possesses an even 
larger polarization, calculated at 100 /iC/cm^ [4l|. Its 
band gap is even lower than BFO, with HSE predict¬ 
ing 1.2 eV In BFO, Bi has oxidation state 3-I-, so 
that its 6s orbital is filled, and the exchange splitting 
of Fe determines the gap. However, in PbNiOs, Pb is 
4-I-, and its 6s-states appear lower in energy than the 
Ni exchange-split bands, resulting in a distinct electronic 
profile. This can be clearly seen in the projected density 
of states (Fig. |4(a)| ): the lowest conduction band is al¬ 
most entirely Pb 6s and O 2p states, while the d-states 
only appear in the valence band and higher in the con¬ 
duction manifold. While this serves to lower the band 
gap dramatically, a further result of this is a Glass coeffi¬ 
cient (Fig. |4(d)| ) over an order of magnitude larger than 
the benchmark value. The current density is modest by 
comparison, though it still exceeds the benchmark, indi¬ 
cating large shift vectors with relatively low absorption. 


HgPbOa and ZnSnOs [6l|, are known to take 
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FIG. 3. The band structures of (a) PbNiOs, |(b)| Mgi/ 2 Zni/ 2 Pb 03 , and (c) LiBiOs. Note the similar, highly dispersive 
conduction band edges. 







FIG. 4. (a) and (b) give the projected density of states for 
PbNiOs. The unfilled half of Cg of the high-spin d® Ni ap¬ 
pears as a sharp peak above the unfilled Pb s-orbitals, which 
have strongly hybridized with oxygen p-orbitals, resulting in 
a low band gap (1.2eV in HSE |44]]1. This material has a 
large [(^ current density response (« 4x benchmark), and a 
very large m Glass coefIicient(« lOx benchmark). 


FIG. 5. Orbital-projected densities of states for 
Mgi/ 2 Zni/ 2 Pb 03 are shown in |(a)| and |(b)| The valence band 
is formed almost entirely from oxygen p-orbitals, and the con¬ 
duction band is hybridized Pb 6s and O 2p-states. This re¬ 
sults in a low band gap (1.2 in HS E), |(c)| high current density 
response (« 4x benchmark), and |(d)| a very large Glass coef¬ 
ficient (« 7x benchmark). 


the ilmenite and LiNbOa structures, respectively. How¬ 
ever, the first is metallic and the second has a high band 
gap and only modest photovoltaic response. We first cal¬ 
culated the response of ZnPbOa, but found it to be bor¬ 
derline metallic, despite promising response; to raise the 
gap, we substituted Mg for half of the Zn. Phonon cal¬ 
culations indicate that the structure is metastable. Once 
again, as seen in Fig. |5(b)[ hybridized Pb 6s states com¬ 
pose the lowest unfilled band. The magnitude of the 
response is quite high, but the current is antiparallel to 
the computed polarization. This is unlike most mate¬ 
rials, including our benchmark materials and the afore¬ 
mentioned PbNiOa, however, we emphasize the ambi¬ 
guity of both polarization (the addition of polarization 


quanta) and structure orientation (designation of A and 
B cations) for these materials. If we compare the two Pb 
compounds with Pb as the B-site in both, not only are 
the structures more similar, but their responses become 
parallel. 

LiBiOs is known to exist in a structure with edge shar¬ 
ing oxygen octahedra [63| . However, our calculations 
place the LiNbOs-type structure - which phonon anal¬ 
ysis reveals to be metastable - nearby in energy, at only 
about O.OI eV per atom higher; additionally, NaBiOs is 
known to take the closely-related ilmenite structure [63| . 
In light of this, we consider it highly possible that the 
LiBiOs can be synthesized in the LiNbOs structure. 

As shown in Fig. 6(a) and Fig. |6(b)[ the electronic 
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(a) (b) 




FIG . 6. The density of states for LiBiOs, shown in |(a)| 
and 1(b) I is dominated by bismuth and oxygen. The band gap 
is set by transitions from O 2p to hybridized Bi 6s states. 
The band gap is modest (1.7eV in HSE). The current density 
response, shown in (c) is quite high (« 8x benchmark), with 
a large [(1 Glass coefficient (« 3x benchmark), indicating 
strong absorption in addition to long shift vectors. 



Direct gap 

X eV 

Max. G 
Xl0"®cm/V 

Max. a 

W/m"^ 

BiFeOs 

2.7 

5 

5 

KBNNO [34| 

1.3(HSE) 

25 

5 

LiAsSe 2 [33] 

1.1 [64] 

98 

30 

PbNiOs 

1.2(HSE) [^ 

50 

25 

Mgi/jZni/jPbOs 

1.2(HSE) 

35 

25 

LiBiOs 

1.7(HSE) 

15 

45 


TABLE II. The band gap and response characteristics of the 
presented materials, along with benchmark BiFeOs and other 
recently proposed bulk photovoltaics (with scissor corrected 
responses) for comparison. 


structure is very similar to the previous two materials. As 
with Pb-containing compounds, the low-lying hybridized 
Bi s states form the lowest unfilled bands, though the 
Bi s proportion is lower than that of Pb s in the afore¬ 
mentioned materials. Possibly as a consequence, the dis¬ 
persion of the conduction band is reduced compared to 
PbNiOa and Mgi/ 2 Zni/ 2 Pb 03 (Fig. |3]), and the BPVE 
response is somewhat different: while the Glass coeffi¬ 
cient is not as large as for the two lead-containing ma¬ 
terials, the photocurrent density is higher, indicating in¬ 
creased absorption. Additionally, the band gap is larger, 
with HSE predicting 1.7-1.8 eV, positioned almost per¬ 
fectly with respect to the visible spectrum for solar en¬ 
ergy conversion. 

It is worth contrasting the response of LiBiOs with 
that of the two lead compounds: the former has a notably 


distinct response, especially near the band edge. This 
can be attributed to the difference in valence band char¬ 
acter; in PbNiOa and Mgi/ 2 Zni/ 2 Pb 03 the valence band 
edge contains considerable density from the secondary 
cation filled d-states. This alters the character of the 
wavefunctions and improves delocalization and response 
magnitude by sharing density, as opposed to lithium’s 
almost completely ionic character. This suggests that in¬ 
clusion of an appropriate dopant with higher electroneg¬ 
ativity may allow for significant tuning of the response 
in LiBiOs. 


We have proposed several polar oxides in the LiNbOs 
structure with strong computed BPVE response and low 
band gaps, summarized in Table [ITl The compositions, 
featuring Pb^+ or Bi^+ cations, were chosen for the ab¬ 
sence of d-states at the band edge. Instead, these mate¬ 
rials have conduction bands formed by low-lying s-states 
hybridized with oxygen p-states. In addition to creating 
significantly lower band gaps, this makes for large, diffuse 
orbitals and strongly delocalized states; combined with 
large polar distortions, they effect significant shift cur¬ 
rent response that is over an order of magnitude higher 
than that previously observed, and roughly double the 
best performing materials previously proposed. Given 
the minimal contributions from the other cations, the 
possibility of tuning the response via composition with¬ 
out altering its fundamental character is strongly sug¬ 
gested. Moreover, in combination with recent demonstra¬ 
tions that careful device construction can dramatically 
improve BPVE performance 13, 23, 3^, these results in¬ 
dicate that BPVE can be much stronger than previously 
thought, bolstering hopes that the phenomenon can be 
successfully exploited. 


S. M. Y. was supported by the Department of En¬ 
ergy Office of Basic Energy Sciences, under grant number 
DE-FG02-07ER46431, and a National Research Coun¬ 
cil Research Associateship Award at the US Naval Re¬ 
search Laboratory. F.Z. was supported by the Office 
of Naval Research under Grant No. N00014-11-1-0664. 
A.M.R was supported by the Office of Naval Research 
under grant N00014-12-1-1033. Computational support 
was provided by the HPCMO of the DoD and NERSC of 
the DOE. 


[1] A. G. Chynoweth, Phys. Rev. 102, 705 (1956). 

[2] F. S. Chen, J. Appl. Phys. 40, 3389 (1969). 

[3] A. M. Glass, D. von der Linde, and T. J. Negran, Appl. 
Phys. Lett. 25, 233 (1974). 

[4] V. M. Fridkin, Crystallog. Rep. 46, 654 (2001). 

[5] R. von Baltz and W. Kraut, Phys. Rev. B 23, 5590 
(1981). 

[6] B. 1. Sturman and V. M. Fridkin, The Photovoltaic and 
Photorefractive Effects in Noncentrosymmetric Materi¬ 
als, edited by G. W. Taylor, Ferroelectricity and Related 



















































6 


Phenomena, Vol. 8 (Gordon and Breach Science Publish¬ 
ers, 1992). 

[7] F. Jermann and J. Otten, J. Opt. Soc. Am. B 10, 2085 
(1993). 

[8] K. Tonooka, P. Poosanaas, and K. Uchino (International 
Society for Optics and Photonics, 1994) pp. 224-232. 

[9] I. Grinberg, D. V. West, M. Torres, G. Gou, D. M. 
Stein, L. Wu, G. Ghen, E. M. Gallo, A. R. Akba- 
shev, P. K. Davies, J. E. Spanier, and A. M. Rappe, 
Nature 503, 509 (2013) 

[10] R. K. Vijayaraghavan, S. C. J. Meskers, M. A. Rahim, 
and S. Das, Ghemical Communications 50, 6530 (2014) 

[11] C. Somma, K. Reimann, C. Flytzanis, T. Elsaesser, and 
M. Woerner, Phys. Rev. Lett. 112, 146602 (2014) 

[12] M. M. Furchi, A. Pospischil, F. Libisch, J. Burgdorfer, 
and T. Mueller, Nano Letters 14, 4785 (2014) 

[13] A. Zenkevich, Y. Matveyev, K. Maksimova, 

R. Gaynutdinov, A. Tolstikhina, and V. Fridkin, 

Phys. Rev. B 90, 161409 (2014) 

[14] F. Zheng, H. Takenaka, F. Wang, N. Z. Koocher, and 
A. M. Rappe, J. Phys. Chem. Lett. 6, 31 (2015) 

[15] M. Ichiki, H. Fume, T. Kobayashi, R. Maeda, 
Y. Morikawa, T. Nakada, and K. Nonaka, Appl. Phys. 
Lett. 87, 222903 (2005). 

[16] M. Qin, K. Yao, and Y. C. Liang, Appl. Phys. Lett. 95, 
022912 (2009). 

[17] Z. J. Yue, K. Zhao, S. Q. Zhao, Z. Q. Lu, X. M. Li, 
H. Ni, and A. J. Wang, J. Phys. D-Appl. Phys. 43, 
015104 (2010). 

[18] G. L. Yuan and J. L. Wang, Appl. Phys. Lett. 95, 252904 
(2009). 

[19] L. Pintilie, V. Stancu, E. Vasile, and 1. Pintilie, J. Appl. 
Phys. 107, 114111 (2010). 

[20] D. W. Gao, H. Zhang, L. A. Fang, W. Dong, F. G. Zheng, 
and M. R. Shen, Appl. Phys. Lett. 97, 102104 (2010). 

[21] S. Y. Yang, J. Seidel, S. J. Byrnes, P. Shafer, C. H. Yang, 
M. D. Rossell, P. Yu, Y. H. Chu, J. F. Scott, J. W. Ager, 
L. W. Martin, and R. Ramesh, Nat. Nanotechnol. 5, 143 
( 2010 ). 

[22] W. Ji, K. Yao, and Y. C. Liang, Adv. Mater. 22, 1763 

( 2010 ). 

[23] M. Alexe and D. Hesse, Nat. Comms. 2, 256 (2011). 

[24] T. Choi, S. Lee, Y. Choi, V. Kiryukhin, and S.-W. 
Cheong, Science 324, 63 (2009). 

[25] R. Nechache, C. Harnagea, S. Licoccia, E. Traversa, 

A. Ruediger, A. Pignolet, and F. Rosei, 

Appl. Phys. Lett. 98, 202902 (2011) 

[26] R. K. Katiyar, P. Misra, F. Mendoza, G. Morell, and 
R. S. Katiyar, Appl. Phys. Lett. 105, 142902 (2014) 

[27] L. Fang, L. You, Y. Zhou, P. Ren, Z. S. Lim, and 
J. Wang, Appl. Phys. Lett. 104, 142903 (2014) 

[28] S. Nakashima, T. Uchida, D. Nakayama, 

H. Fujisawa, M. Kobune, and M. Shimizu, 
Jpn. J. Appl. Phys. 53, 09PA16 (2014) 

[29] S. Gupta, R. Medwal, T. B. Limbu, R. K. Katiyar, S. P. 
Pavunny, M. Tomar, G. Morell, V. Gupta, and R. S. 
Katiyar, Appl. Phys. Lett. 107, 062902 (2015) 

[30] R. L. Gao, C. L. Fu, W. Cai, G. Chen, X. L. 
Deng, H. W. Yang, J. R. Sun, and B. G. Shen, 
Materials Science Forum 815, 176 (2015) 

[31] R. Nechache, C. Harnagea, S. Li, L. Carde¬ 
nas, W. Huang, J. Chakrabartty, and F. Rosei, 
Nature Photonics 9, 61 (2014) 

[32] L. Jiang, 1. Grinberg, F. Wang, S. M. 


Young, P. K. Davies, and A. M. Rappe, 
Phys. Rev. B 90, 075153 (2014) 

[33] J. A. Brehm, S. M. Young, F. Zheng, and A. M. Rappe, 
The Journal of Chemical Physics 141, 204704 (2014) 

[34] F. Wang and A. M. Rappe, 
Phys. Rev. B 91, 165124 (2015) 

[35] J. E. Sipe and A. 1. Shkrebtii, Phys. Rev. B 61, 5337 

( 2000 ). 

[36] A. Bhatnagar, A. R. Chaudhuri, Y. H. Kim, D. Hesse, 
and M. Alexe, Nature Communications 4, 2835 (2013) 

[37] S. M. Young and A. M. Rappe, 
Phys. Rev. Lett. 109, 116601 (2012) 

[38] A. M. Rappe, K. M. Rabe, E. Kaxiras, and J. D. 
Joannopoulos, Phys. Rev. B Rapid Comm. 41, 1227 
(1990). 

[39] N. J. Ramer and A. M. Rappe, Phys. Rev. B 59, 12471 
(1999). 

[40] S. M. Young, F. Zheng, and A. M. Rappe, 
Phys. Rev. Lett. 109, 236601 (2012) 

[41] X. F. Hao, A. Stroppa, S. Picozzi, A. Filippetti, and 
C. Franchini, Phys. Rev. B 86, 014116 (2012) 

[42] P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, 
C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, M. Cococ- 
cioni, 1. Dabo, A. D. Corso, S. de Gironcoli, S. Fabris, 
G. Fratesi, R. Gebauer, U. Gerstmann, C. Gougoussis, 
A. Kokalj, M. Lazzeri, L. Martin-Samos, N. Marzari, 
F. Mauri, R. Mazzarello, S. Paolini, A. Pasquarello, 
L. Paulatto, C. Sbraccia, S. Scandolo, G. Sclauzero, A. P. 
Seitsonen, A. Smogunov, P. Umari, and R. M. Wentz- 
covitch, J. Phys.: Condens. Matter 21, 395502 (2009). 

[43] J. Heyd, G. E. Scuseria, and M. Ernzerhof, The Journal 
of Chemical Physics 118, 8207 (2003). 

[44] Z. H. Levine and D. C. Allan, 
Phys. Rev. Lett. 63, 1719 (1989) 

[45] F. Nastos, B. Olejnik, K. Schwarz, and J. E. Sipe, Phys. 
Rev. B 72, 045223 (2005). 

[46] K. A. Rabe, M. Dawber, C. Lichtensteiger, C. H. Ahn, 
and J. M. Triscone, Topics in Applied Physics 105, 1 
(2007). 

[47] J. B. Neaton, C. Ederer, U. V. Waghmare, N. A. Spaldin, 
and K. M. Rabe, Phys. Rev. B 71, 014113 (2005) 

[48] S. V. Levchenko and A. M. Rappe, 

Phys. Rev. Lett. 100, 256101 (2008) 

[49] C. Baeumer, D. Saldana-Greco, J. M. P. Martirez, 

A. M. Rappe, M. Shim, and L. W. Martin, 

Nature Communications 6, 6136 (2015) 

[50] Y. Inoue, Energy Environ. Sci. 2, 364 (2009) 

[51] L. G. Reznik, A. A. Anikiev, B. S. Umarov, and J. F. 
Scott, Ferroelectrics 64, 215 (1985) 

[52] A. Anikiev, L. G. Reznik, B. S. Umarov, and J. F. Scott, 
Ferroelectr. Lett. 3, 89 (1985) 

[53] A. Dhar and A. Mansingh, 

J. Appl. Phys. 68, 5804 (1990) 

[54] G. Catalan and J. F. Scott, Adv. Mater. 21, 2463 (2009). 

[55] J. F. Ihlefeld, N. J. Podraza, Z. K. Liu, R. C. Rai, X. Xu, 
T. Heeg, Y. B. Chen, J. Li, R. W. Collins, J. L. Musfeldt, 
X. Q. Pan, J. Schubert, R. Ramesh, and D. G. Schlom, 
Appl. Phys. Lett. 92, 142908 (2008). 

[56] J. Seidel, D. Y. Fu, S. Y. Yang, E. Alarcon-Llado, J. Q. 
Wu, R. Ramesh, and J. W. Ager, Phys. Rev. Lett. 107, 
126805 (2011). 

[57] W. Ji, K. Yao, and Y. C. Liang, Phys. Rev. B 84, 094115 

( 2011 ). 

[58] Y. Inaguma, K. Tanaka, T. Tsuchiya, D. Mori, T. Kat- 





7 


sumata, T. Ohba, K.-i. Hiraki, T. Takahashi, and 
H. Saitoh, J. Am. Chem. Soc. 133, 16920 (2011) 

[59] X. F. Hao, A. Stroppa, P. Barone, A. Fil- 

ippetti, C. Franchini, and S. Picozzi, 

New Journal of Physics 16, 015030 (2014) 

[60] A. Sleight and C. Prewitt, 

J. Solid State Chem. 6, 509 (1973) 

[61] Y. Inaguma, M. Yoshida, and T. Katsumata, 
J. Am. Chem. Soc. 130, 6704 (2008) 

[62] H. Gou, F. Gao, and J. Zhang, 


Comput. Mater. Sci. 49, 552 (2010) 

[63] T. Takei, R. Haramoto, Q. Dong, N. Ku- 
mada, Y. Yonesaki, N. Kinomura, T. Mano, 
S. Nishimoto, Y. Kameshima, and M. Miyake, 
J. Solid State Chem. 184, 2017 (2011) 

[64] T. K. Bera, J. I. Jang, J.-H. Song, C. D. Malliakas, 
A. J. Freeman, J. B. Ketterson, and M. G. Kanatzidis, 
J. Am. Chem. Soc. 132, 3484 (2010) 



